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ABSTRACT 

A  one  dimensional  (1-D),  isothermal  model  for  a  direct  methanol  fuel  cell 
(DMFC)  is  presented.  This  model  accounts  for  the  kinetics  of  the  multi-step  methanol 
oxidation  reaction  at  the  anode.  Diffusion  and  crossover  of  methanol  are  modeled  and  the 

mixed  potential  of  the  oxygen  cathode  due  to  methanol  crossover  is  included.  Kinetic  and 

2 

diffusional  parameters  are  estimated  by  comparing  the  model  to  data  from  a  25  cm 
DMFC.  This  semi-analytical  model  can  be  solved  rapidly  so  that  it  is  suitable  for 
inclusion  in  real-time  system  level  DMFC  simulations. 

Keywords:  DMFC;  Analytical  Model;  Mixed  Potential;  Methanol  Crossover 
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INTRODUCTION 

Direct  Methanol  Fuel  Cells  (DMFCs)  are  currently  being  investigated  as 
alternative  power  source  to  batteries  for  portable  applications  because  they  can  offer 
higher  energy  densities.  However,  two  factors  limit  the  performance  of  DMFC  systems: 
crossover  of  methanol  from  anode  to  cathode  and  the  slow  kinetics  of  the  electrochemical 
oxidation  of  methanol  at  the  anode. 

The  crossover  of  methanol  lowers  the  system  efficiency  and  decreases  cell 
potential  due  to  corrosion  at  the  cathode.  Figure  1  illustrates  the  electrochemistry  and 
transport  phenomena  in  DMFCs.  Electrochemical  oxidation  of  methanol  occurs  at  both 
anode  and  cathode,  but  corrosion  current  at  the  cathode  produces  no  usable  work.  Several 
experimental  and  modeling  studies  have  characterized  methanol  crossover  in  DMFCs  [1- 

4]- 

The  kinetics  of  DMFCs  are  complicated  because  the  reaction  mechanism  involves 
adsorption  of  methanol  and  several  reaction  steps  including  the  oxidation  of  CO.  Figure  2 
shows  a  possible  network  of  reaction  pathways  by  which  the  electrochemical  oxidation  of 
methanol  occurs.  Catalysis  studies  have  attempted  to  analyze  possible  reaction  pathways 
to  find  the  main  pathway  of  methanol  oxidation  [6-8].  Most  studies  conclude  that  the 
reaction  can  proceed  according  to  multiple  mechanisms.  However,  it  is  widely  accepted 
that  the  most  significant  reactions  are  the  adsorption  of  methanol  and  the  oxidation  of 
CO.  Follows  is  a  simplified  reaction  mechanism  that  will  be  used  in  this  paper  to  model 
performance  of  DMFCs. 

Cl  IX)  1 1  +  Site  ( Cl  1 ,01 1  )ads  ( 1 ) 


(CHfiH^^CO^+W+Ae 


(2) 
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(CO)ads  +  H20  ->  C02  +  2H+  +  2e  (3) 

This  mechanism  is  similar  to  the  mechanism  used  by  Meyers  and  Newman  [9],  but  does 
not  segregate  the  electrochemical  oxidation  of  water  reaction  from  the  electrochemical 
oxidation  of  CO.  This  assumption  does  not  change  the  kinetic  expression  appreciably  and 
is  applicable  for  Pt-Ru  catalysts  where  the  oxidation  of  water  on  Ru  occurs  much  faster 
than  the  oxidation  of  CO. 

The  model  presented  in  this  paper  seeks  to  provide  a  one  dimensional  (1-D), 
isothermal  model  of  a  DMFC  that  allows  rapid  prediction  of  polarization  data  and  gives 
insight  into  mass  transport  phenomena  occurring  in  the  cell.  Models  currently  in  the 
literature  leave  out  effects  important  for  predicting  full  cell  performance  or  include 
physical  detail  that  encumbers  the  model  and  complicates  its  solution.  Baxter  et  al.  [10] 
developed  a  model  for  the  DMFC  anode  which  considers  diffusion  of  C02,  FFO  and 
methanol  in  the  anode,  but  neglects  the  effects  of  the  cathode  and  thus  does  not  capture 
the  effects  of  methanol  crossover.  They  also  used  Butler- Volmer  kinetics  to  describe  the 
electrochemical  oxidation  of  methanol.  Meyers  and  Newman  [9]  develop  a  kinetic 
expression  similar  to  the  one  used  in  this  paper  and  provide  parameters  for  the  cathode 
reaction.  The  level  of  detail  included  in  their  membrane  model  and  the  resulting  non¬ 
linearity  of  their  equations  make  the  solution  of  their  model  difficult.  Kulikovsky  [11] 
solved  an  analytical  model  for  the  fuel  cell  anode  to  predict  the  anodic  overpotential. 
However,  the  model  presented  could  only  be  solved  in  the  limits  of  low  current  or  high 
current  and  did  not  allow  prediction  of  full  cell  polarization  behavior.  Wang  and  Wang 
[12]  used  a  Computational  Fluid  Dynamics  (CFD)  model  to  investigate  a  full  DMFC  fuel 
cell.  This  analysis  included  two-phase  flow  effects  in  the  backing  layers  (BLs)  but  used  a 
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non-intuitive  transition  in  kinetics  at  a  certain  concentration  to  describe  the  complex 
methanol  oxidation  reaction.  Norlund  and  Lindbergh  [13]  develop  an  anode  model  that 
neglects  the  effects  of  methanol  crossover  and  the  cathode.  Their  model  also  assumes  a 
flooded  agglomerate  model  of  the  anode  catalyst  layer  (ACL)  that  assumes  a  specific 
geometry  for  all  reaction  sites. 

EXPERIMENTAL 

Cell  Preparation 

Tests  were  performed  on  a  25  cm2  fuel  cell  from  Fuel  Cell  Technologies.  The 
membrane  electrode  assembly  (MEA)  was  constructed  from  a  Nafion®  117  membrane 
with  E-TEK  40%  Platinum/C  gas  diffusion  electrodes  prepared  according  to  the  decal 
method  of  Wilson  [14].  The  anode  loading  was  3  mg/cm  of  1:1  Pt/Ru  catalyst  and  the 
cathode  loading  was  1  mg/cnL  of  Pt.  Tests  were  conducted  using  an  890C  load  cell  from 
Scribner  Associates  Inc.  with  a  methanol  fuel  system.  The  cell  was  broken  in  by  running 
for  3  hours  under  a  5  A  load  with  a  40  mL/min  flow  of  1  M  feed  to  the  anode  and  50 
mL/min  flow  of  dry  oxygen  on  the  cathode.  The  cell  temperature  and  inlet  temperatures 
were  70°C.  All  reagents  were  certified  as  ultra  high  purity. 

Testing 

Prior  to  running  tests  with  a  given  concentration  of  methanol  the  system  was 
flushed  with  1.5  L  of  methanol  solution.  The  flow  rates  for  the  anode  and  cathode  were 
then  set  to  those  necessary  to  maintain  10/5  stoichiometric  excess  ratios  on  the  anode  and 
cathode.  The  minimum  flow  rate  for  all  experiments  was  10  mL/min  on  the  anode  and  50 
mL/min  on  the  cathode.  The  cell  was  next  run  under  a  load  of  0.25 A  for  10  minutes  or 


until  the  voltage  reached  steady-state.  The  load  was  set  to  0  A  for  10  minutes  or  until  the 
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voltage  arrived  at  its  steady  open  circuit  value.  Polarization  curves  were  run  in  current 
scan  mode  with  0.04  A/point  and  150  seconds/point. 

MODEL  DEVELOPMENT 

Assumptions 

The  model  presented  here  accounts  for  concentration  variations  of  methanol 
across  the  anode  backing  layer  (ABL),  anode  catalyst  layer  (ACL),  and  membrane. 
Figure  3  presents  a  schematic  diagram  of  the  layers  considered  in  the  model  illustrating 
several  assumptions.  The  assumptions  used  in  this  model  are 

1 .  Steady-state. 

2.  Variations  in  only  one  spatial  Cartesian  coordinate  (i.e.,  across  the  MEA). 

3.  Convective  transport  of  methanol  is  negligible. 

4.  Pressure  gradient  across  the  layers  are  negligible. 

5.  Isothennal  conditions. 

6.  All  physical  properties  are  considered  constant. 

7.  Only  liquid  phase  is  considered.  This  means  that  carbon  dioxide  remains 
dissolved  in  solution. 

8.  Solutions  are  considered  ideal  and  diluted. 

9.  Local  equilibrium  at  interfaces  can  be  described  by  a  partition  function. 

10.  The  ACL  is  assumed  to  be  a  macro-homogeneous  porous  electrode  and  thus  the 
reaction  in  this  layer  is  modeled  as  a  homogeneous  reaction. 

11.  Anode  kinetics  can  be  described  by  the  step  mechanism,  Eq.  (1)  -  (3),  with  a  rate 
expression  similar  to  the  one  obtained  by  Meyers  and  Newman  [9]. 

12.  The  anodic  overpotential  is  constant  throughout  the  ACL. 
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13.  Cathode  kinetics  can  be  described  by  Tafel  expression  with  no  mass  transfer 
limitations. 

Applying  these  assumptions,  the  mass  transport  equations  are  developed  and 
combined  with  the  kinetic  equations  in  order  to  calculate  the  cell  voltage,  which  can  be 
expressed  as: 


VCen  =  u°2  -  UMeOH  -r]c-rjA 


^ M I  Cell 

K 


(4) 


where  U°2  and  U MeOH  are  the  thermodynamic  equilibrium  potential  of  oxygen  reduction 
and  methanol  oxidation  respectively,  ij(  and  //,  are  the  cathode  and  anode 

overpotentials,  respectively,  and  the  last  tenn  in  Eq.  (4)  represents  the  ohmic  drop  across 
the  membrane. 

Governing  Equations  and  Boundary  Conditions-Anode 

The  anode  overpotential  is  obtained  by  first  obtaining  the  concentration  profiles 
across  the  various  regions  of  the  MEA. 

Anode  Backing  Layer 

The  differential  mass  balance  for  methanol  in  the  ABL  is 


dNB 

ulyMeOH,z  _q 

dz 


(5) 


Assuming  Fickian  diffusion  [15]  of  methanol  with  an  effective  diffusivity  DB  in  the  ABL 
phase,  the  methanol  flux  gives 


N, 


MeOH,z  B 


r  B 

_£)  aCMeOH 


dz 


(6) 


Substitution  of  Eq.  (6)  into  Eq.  (5)  gives  the  governing  equation  for  methanol  in  the  ABL 


as 
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i2  B 

'  CMeOH  _  Q 

J_2 


The  boundary  conditions  for  Eq.  (7)  are  illustrated  in  Fig.  3.  It  is  assumed  that 
concentration  at  the  flow-channel/ ARE  interface  is  given  by  the  bulk  concentration  in  the 
flow  channel.  The  concentration  at  the  ABL/ACL  interface  is  given  by  assuming  local 
equilibrium  with  a  partition  coefficient  Kl . 


Atz  =  0:  cMeOH=cb 


A  J.  .  B  B  TS-  A 

At  z  —  Zj .  cMeOH  —  Cj  —  KIcI 


Membrane 


The  differential  mass  balance  for  methanol  in  the  membrane  is 


dNM 

MeOH,z  _  q 


The  transport  of  methanol  in  the  membrane  is  governed  by  diffusion  and  electro-osmotic 


drag.  The  flux  equation  can  be  written  as 


*  MeOH,z 


i  M 

...  _r)  aCMeOH  ,  e 
‘-'M  ,  "l"  bA; 

dz 


where  Du  and  gMe()ri  are  the  effective  diffusion  and  the  electro-osmotic  drag  coefficients 

of  methanol  respectively.  The  electro-osmotic  drag  coefficient  is  defined  as  the  number 
of  methanol  molecules  dragged  by  a  hydrogen  ion  moving  in  the  membrane.  Substitution 
of  Eq.  (11)  into  Eq.  (10)  gives  the  governing  equation  for  methanol  in  the  membrane  as 


'-'MeOH  _  0 

dz 2 


The  boundary  conditions  for  Eq.  (12)  are  illustrated  in  Fig.  3.  It  is  assumed  that  all 
the  methanol  crossing  the  membrane  reacts  at  the  cathode  creating  a  very  low 
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concentration  at  the  membrane/cathode-layer  interface.  The  concentration  at  the 
ACL/membrane  interface  is  given  by  assuming  local  equilibrium  with  a  partition 
coefficient  Ku . 


At  z  =  zu  :  cZOH=c?=Kllc£  (13) 

At  z  =  zm  :  c^eOH  ~  0  (14) 

Anode  Catalyst  Layer 

The  methanol  oxidation  reaction  at  the  anode  is  considered  homogeneous.  The 
differential  mass  balance  for  methanol  in  the  ACL  is 


dNA 

UJV  MeOH,z 


dz 


rMeOH 
M  MeOH 


(15) 


where  the  molar  consumption  rate  [rMe0H  IMUeOH)  is  related  to  the  volumetric  current 
density  j  as 


rMeOH 
M MeOH 


6  F 


(16) 


The  current  density  expression  for  methanol  oxidation  is  taken  from  Meyers  and 
Newman  [9]  as 


j  =  aIi 


MeOH 

0,ref 


kc 


^ MeOH 


aArjAF 

MeOH _ e  RT 

aAlAF 

+  Ae  RT 


(17) 


where  a  is  the  specific  surface  area  of  the  anode,  I^H  is  the  exchange  current  density, 
and  k  and  X  are  constants. 

The  methanol  flux  in  the  ACL  with  an  effective  diffusivity  D4  is  given  by  a 
similar  expression  as  showed  for  the  ABL. 
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(18) 

dz 

Substitution  of  Eq.  (16)  and  (18)  into  Eq.  (15)  gives  the  governing  equation  for  methanol 
in  the  ACL  as 

Z)^42»=i  (19) 

dz  6  F 

The  boundary  conditions  for  Eq.  (19)  are  illustrated  in  Fig.  3.  The  methanol 
concentration  at  the  interfaces  is  given  as 


At  z  =  z1\ 

§ 

II 

Ci 

(20) 

At  z  = zu  : 

A  _  A 

CMeOH  ~  CII 

(21) 

The  concentrations  given  in  Eq.  (20)  and  (21)  are  related  to  the  concentrations  at  the  ABL 
and  the  membrane  through  Eq.  (9)  and  (13).  These  concentrations  can  be  determined 
from  jump  mass  balances  [15]  at  the  ABL/ACL  and  ACL/membrane  interfaces,  yielding 

At  z  =  Zj\  NBMeOHyZ  =  NAMeOH'Z  (22) 


At  z  =  zu  :  Na  =  N™ 


Analytical  Solution-Anode 

The  solution  to  Eq.  (7)  -  (9)  is 


B  _  K ]C f  Cb 

CMeOH  ~  e  Z  +Cb 


The  solution  to  Eq.  (12)  -  (14)  is 


cM  =  K  CA\^B+^A _ -  +  \ 

^ MeOH  Il^ll  c  T  1 

k  0*4 


The  solution  to  Eq.  (19)  -  (21)  is: 
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‘Cell  _2 


MeOH 


12  FSaDa 


z1+Clz  +  C2 


(26) 


where 


C, 


_  4  -cf  fe„44+4) 

*a  UF5aDa 


(27) 


and 


C2=cf- 


(4~cf)SB  |  ICellSB(SB+SA) 


S, 


12  F8aDa 


(28) 


From  the  solutions  above  the  fluxes  in  each  phase  can  be  obtained  via  Eq.  (6),  (11),  and 
(18).  The  fluxes  are  then  evaluated  at  the  respective  interfaces  to  obtain  two  expressions 
in  terms  of  cf  and  cf  from  Eq.  (22)  and  (23).  One  may  ultimately  show  that 


sadmku 


f  I  S'' 

/)  „  „  1CelluB 


V 


12F 


+  8mDa 


(  IS'' 

D  C  -0  +  6(5  yCellUB 


DbK ,  (, 8aDmKu  +  8mDa  )  +  SBDADMKh 


(29) 


4  = 


L 


Ir 


DADBcb  -  SaDbKj  (1  + 12^ )  -f'L  _  SbDa  (1  +  6{MeOH )  ^ 
v  2 nF  6F  j 

DbK2  ( 8aDmKu  +  SmDa  )  +  SbDaDmKu 


(30) 


Finally,  the  concentration  profile  given  by  Eq.  (26)  is  substituted  into  the  kinetic 
expression,  Eq.  (17),  integrated,  and  equated  to  the  cell  current  giving 


SB+$A 


I  =  f  aIM 

1  Cell  J  U1 0, 


MeOH 

ref 


krA 

'vt' MeOH 


<* a1af 

c A  +  2.P  RT 

^  MeOH  ^ 


aArliP 

?  RT  dz 


(31) 


Assuming  rjA  is  constant  (assumption  12),  Eq.  (31)  is  used  to  obtain  rjA  for  a  given  value 


of  feu  ■ 


Cathode 


Tafel  kinetics  with  first  order  oxygen  concentration  dependence  is  employed  to 
describe  the  oxygen  reduction  at  the  cathode. 
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c  ac1cF 
C02,ref 


(32) 


where  Ileak  is  the  leakage  current  density  due  to  the  oxidation  of  methanol  crossing  the 
membrane.  The  leakage  current  density  can  be  written  as 


/  =6 FNm 

1  leak  iy  MeOH,z 


(33) 


where  N^eOH  z  is  obtained  from  Eq.  (11).  Equation  (32)  is  then  used  to  obtain  r/c  for  a 
given  value  of  ICell . 

Finally,  the  anode  and  cathode  overpotentials  are  substituted  into  Eq.  (4)  to  give 
give  VCell  for  a  given  value  of  ICelI . 


RESULTS  AND  DISCUSSION 


Experimental  and  modeling  results  of  polarization  behavior  for  0.05M,  0.1M, 
0.2M,  and  0.5M  methanol  solutions  are  shown  in  Fig.  4.  The  limiting  current  densities 
predicted  by  the  model  are  very  close  to  experimental  values.  The  model  predictions  for 
conditions  near  open  circuit  voltage  show  the  largest  errors  with  experimental  values. 
This  disagreement  could  be  due  to  the  fact  that  concentration  and  temperature  effects  on 
the  thermodynamic  potentials  of  the  electrodes  were  neglected.  Methanol  polarization 
data  above  0.5M  could  not  be  modeled  with  the  same  set  of  kinetic  and  transport 
parameters  as  was  used  for  the  cases  shown  in  Fig.  4.  Trends  in  the  predicted  and 
modeled  polarization  curves  in  Fig.  4  are  similar  to  those  shown  for  0.2M  and  0.5M  in 
Wang  and  Wang  [12].  However,  the  limiting  current  densities  Wang  and  Wang  [12] 
predict  are  higher  than  those  in  Fig.  4.  In  their  paper,  they  contend  that  high  current 
densities  in  DMFCs  can  be  explained  by  the  possibility  of  gas  phase  transport. 
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The  modeling  parameters  used  are  listed  in  Table  1.  Transport  parameters  agree 
well  with  literature  values.  The  specific  area  (a)  and  the  anode  and  cathode  transfer 
coefficients  can  change  due  to  electrode  properties  and  were  adjusted  to  fit  the  model  to 
the  experimental  data.  It  was  found  that  around  the  parameter  set  listed  in  Table  1  certain 
parameters  could  be  adjusted  simultaneously  and  the  resulting  fit  did  not  alter  the 
polarization  curves  significantly.  One  example  is  that  increasing  the  exchange  current 
density  while  increasing  X  produced  nearly  equivalent  curves.  For  this  reason,  all 
parameters  in  Table  1  are  listed  only  to  two  significant  digits.  For  the  model  development 
the  methanol  electro-osmotic  drag  coefficient  was  assumed  to  be  a  constant  value,  but 
when  solving  the  model  the  methanol  electro-osmotic  drag  coefficients  was  estimated  at 
every  point  on  the  polarization  curve  according  to  the  equation  in  Table  1. 

Figure  5  shows  predicted  concentration  profiles  across  the  anode  and  membrane 
for  the  four  concentrations  at  15  mA/cm".  At  this  condition  a  cell  operating  with  a  0.05M 
bulk  methanol  concentration  is  in  the  mass  transfer  limited  region  while  the  0.1M,  0.2M, 
and  0.5M  concentrations  are  in  the  region  limited  by  the  oxidation  of  CO  on  the  catalyst 
surface.  The  concentration  for  the  0.05M  case  in  the  catalyst  layer  is  very  low  at  this 
current  density  similar  to  what  should  be  expected.  The  concentration  profile  across  the 
catalyst  layer  appears  to  be  nearly  constant  for  the  0.05M,  0.1M,  and  0.2M 
concentrations.  The  0.5M  concentration  has  a  larger  drop  in  concentration  across  the 
catalyst  layer  due  to  a  higher  rate  of  methanol  crossover,  but  the  value  is  still  relatively 
constant.  The  assumption  that  the  methanol  concentration  in  the  ACL  is  constant  is  most 
valid  close  to  the  limiting  current  density  where  the  methanol  concentration  is  the  lowest, 
thus  reducing  the  amount  of  methanol  crossover. 
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Figure  6  shows  calculations  of  the  methanol  crossover  predicted  by  the  model  as  a 
function  of  current  density.  At  the  cathode  the  methanol  that  crosses  the  membrane  is 
oxidized  in  a  corrosion  reaction.  The  leakage  current  cannot  be  used  to  do  work. 
Expressing  the  methanol  crossover,  as  in  Fig.  6,  in  terms  of  the  leakage  current  gives  a 
more  tangible  understanding  of  the  loss  in  efficiency  due  to  methanol  crossover.  The 
leakage  current  can  be  reduced  by  running  the  cell  at  low  methanol  concentrations  and 
high  current  densities.  Thus  to  reduce  crossover  running  at  lower  concentrations  of 
methanol  may  be  advantageous.  The  leakage  currents  calculated  in  this  paper  are  similar 
to  those  calculated  by  Wang  and  Wang  [12].  It  should  be  noted  that  the  leakage  current 
goes  to  zero  at  the  limiting  current  value  for  all  concentrations.  This  provides  a  check  that 
our  transport  equations  are  giving  a  physically  meaningful  concentration  profile. 

CONCLUSIONS 

A  semi-  analytical,  1-D,  isothennal  model  of  a  DMFC  has  been  developed.  Using 
reasonable  transport  and  kinetic  parameters  the  model  fits  well  to  experimental 
polarization  data.  The  model  allows  prediction  of  concentration  profiles  in  the  anode  and 
membrane  as  well  as  estimating  methanol  crossover.  The  solution  time  is  less  than  1 
minute. 
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NOMENCLATURE 

specific  surface  area  of  the  anode,  cm'1 

3 

bulk  concentration  of  methanol  in  the  flow  channel,  mol/cm 

3 

concentration  of  methanol  at  the  ABL/ACL  interface,  mol/cm 

3 

concentration  of  methanol  at  the  ACL/membrane  interface,  mol/cm 

-5 

concentration  of  methanol,  mol/cm 
concentration  of  oxygen,  mol/cm 

3 

total  concentration  in  the  ABL,  mol/cm 

2 

effective  diffusion  coefficient  of  methanol  in  the  ACL,  cm"/s 
effective  diffusion  coefficient  of  methanol  in  the  ABL,  cm"/s 
effective  diffusion  coefficient  of  methanol  in  the  membrane,  cm“/s 

Faraday’s  constant,  96,487  C/equiv 
cell  current  density,  A/cm- 

2 

leakage  current  density  due  to  methanol  crossover,  A/cm 
exchange  current  density  of  methanol,  A/cm2 
exchange  current  density  of  oxygen,  A/cm- 

3 

volumetric  current  density,  A/cm 
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k  constant  in  the  rate  expression  (Eq.  21),  dimensionless 

MMeOH  molecular  weight  of  methanol,  g/mol 

N.  MeOH  z  component  of  methanol  molar  flux,  mol/(cm2  s) 

R  gas  constant,  8.314  J/(mol  K) 

•5 

rMeOH  rate  of  consumption  of  methanol  by  homogeneous  reaction,  g/  (cm  s) 

T  temperature,  K 

UMeOH  thennodynamic  equilibrium  potential  of  methanol  oxidation,  V 

U°2  thennodynamic  equilibrium  potential  of  oxygen  oxidation,  V 

VCell  cell  voltage,  V 

xMeOH  mole  fraction  of  methanol,  mol/mol 

z  coordinate  direction  nonnal  to  the  anode,  cm 

Greek 

aA  anodic  transfer  coefficient 

ac  cathodic  transfer  coefficient 

SA  ACL  thickness,  cm 

SB  ABL  thickness,  cm 

SM  membrane  thickness,  cm 

1a 


anode  overpotential,  V 
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rjc  mix  overpotential  at  the  cathode,  V 

k  ionic  conductivity  of  the  membrane,  S/cm 

3 

A  constant  in  the  rate  expression  (Eq.  21),  mol/cm 

‘zMeOH  electro-osmotic  drag  coefficient  of  methanol 

Subscripts 
A  ACL 

B  ABL 

b  bulk 

Cell  cell 

I  ABL/ACL  interface 

II  ACL/membrane  interface 

III  membrane/cathode  layer  interface 

M  membrane 

MeOH  methanol 

O?  oxygen 

z  z-direction 

Superscripts 
A  ACL 


B 


ABL 
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M 

membrane 

MeOH 

methanol 

o2 

oxygen 
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Cathode:  -0,  +6 H*  +6e  ->3 A/,0 
2  - 

Cathode 

corrosion  CH}OH  +  //,0  — »C0,+6//  +6e 
reaction: 


Overall: 


1 0,  +  C//,0//  — >  C0,  +  2  A/,0 


Figure  1.  Schematic  of  a  DMFC. 
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Figure  2.  Reaction  pathways  of  methanol  oxidation  [5]. 
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Figure  3.  Schematic  of  the  DMFC  layers  considered  in  the  model. 


0.05  M  0.1 


0 


0  0.02 


Figure  4 


Methanol  concentration  (mol/L) 


Brenda  L.  Garcia  26 


Leakage  current  (A/cni 


Brenda  L.  Garcia  27 
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Figure  6.  Methanol  crossover  for  different  methanol  bulk  concentrations. 
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Table  1.  Parameter  Values. 


Parameter 

Value 

Ref. 

a 

1000  cm2 

Assumed 

da 

2436f  — 1 — — )  cm2 

2.8x10  1353  t] 

s 

Scott  et  al.  [16] 

Db 

8.7x1  O'6  cm2/s 

Assumed 

Dm 

4.9xl0-6e2436fe"^  Cnr 

s 

Scott  et  al.  [16] 

jMeOH 

1 0,ref 

35570f  1  O 

9.425x1 0~3e  *  ^353  T>  , 

cm 

Wang  and  Wang  [12] 

Itf 

73200/'  1  lX 

4.222x1 0~3e  R  ^353  T>  , 

cm 

Parthasarathy  et  al.  [17] 

KI 

0.8 

Baxter  et  al.  [10] 

Kn 

0.8 

Baxter  et  al.  [10] 

k 

7.5xl0‘4 

Assumed 

T 

343.15  K 

jjMeOH 

0.03  V 

Wang  and  Wang  [12] 

U° 2 

1.24  V 

Wang  and  Wang  [12] 

°A 

0.52 

Assumed 

°C 

1.55 

Assumed 

8a 

0.0023  cm 

8b 

0.015  cm 
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5m 

0.018  cm 

K 

0.036  S/cm 

Assumed 

X 

2.8xl0'9  mol/cm3 

Assumed 

^ MeOH 

2  5x 

MeOH 

Ren  et  al.  [18] 

